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ABSTRACT 

A new type of gravitational microlensing experiment toward a field where 
stars are not resolved is being developed observationally and theoretically: 
pixel lensing. When the experiment is carried out toward the M31 bulge area, 
events may be produced both by Massive Compact Halo Objects (MACHOs) 
in our Galactic halo and by lenses in M31. We estimate that ~ 10 — 15% of the 
total events are caused by Galactic halo MACHOs assuming an all-MACHO 
halo. If these Galactic events could be identified, they would provide us with an 
important constraint on the shape of the halo. We test various observables that 
can be used for the separation of Galactic halo/M31 events. These observables 
include the Einstein time scale, the effective duration of an event, and the 
flux at the maximum amplification, but they cannot be used to separate each 
population events. However, we find that most high maximum-flux Galactic 
halo events can be isolated through a satellite-based measurement of the flux 
difference caused by the parallax effect. For the detection of the flux difference, 
it is required to monitor events with an exposure time of ~ 20 min by a 0.5 
m telescope mounted on a satellite. Such observations could be carried out 
as a minor component of a mission aimed primarily at events seen toward 
the Galactic bulge and Large Magellanic Cloud. In addition, proper motion 
can be used to isolate Galactic halo/M31 events, but only for ~ 5% of high 
signal-to-noise ratio M31 events and only 1% of Galactic halo events. 
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1. Introduction 

Crotts (1992) and Baillon et al. (1993) have begun a pioneering search for lensing 
events of unresolved stars in M31 by monitoring individual pixels for the time-dependent 
flux induced by one of the many unresolved stars they contain: 'pixel lensing' for short. 
For a classical lensing event, the light curve is obtained by subtracting the reference flux, 
B, from the amplified flux, F: 

F 0:i (A-l) = F-B; B = F ,i + B rcs , (1.1) 

where Fq^ is the flux of the lensed star before or after the event and B rcs = J2j^i Foj is the 
residual flux from stars that are not participating in the event. Here the amplification is 

A(x) = x(x f + + 4 2 )1/2 , x = + u\t - i )T /2 , (1.2) 

where u> = i" 1 is the inverse Einstein ring crossing time, j3 is the dimensionless impact 
parameter, and to is the time at the maximum amplification. In an uncrowded field where 
the blending of stars is not important (B ms <C Fq^), one can directly measure the values 
of F and B = Fq^. On the other hand, it is impossible to measure the absolute values of 
F and B for the case of pixel lensing since stars cannot be resolved. However, the light 
curve, Fq,i{A — 1), can still be obtained by subtracting the reference image from those in 
which a lensing event is in progress (Ciardullo et al. 1990; Tomaney & Crotts 1996). 

Pixel lensing has several advantages over classical lensing. First of all, one is not 
restricted to observing fields where stars are resolvable. Second, one can detect O(10 2 ) 
events from a season- long observation of M31 in a single field with moderate resources 
because a large number of stars are monitored (Han 1996). Finally, one looks through 
the Galactic halo when an external galaxy is monitored for pixel lensing; some fraction 
of events may be caused by the Galactic Massive Compact Halo Objects (MACHOs). 
If the Galactic MACHOs can be distinguished from events caused by the lenses in the 
external galaxy, pixel lensing will considerably increase the number of Galactic MACHO 
events. In addition to the increase in the number of events, one can use the data toward 
these additional lines of sight to constrain the shape of the Galactic halo. However, events 
caused by one population are just contaminants in the study of the other population unless 
they are separated. The problem becomes more serious when the fraction of contaminating 
events is large. It is therefore important to determine the fraction of events caused by each 
population and to develop methods for distinguishing them from one another. 

We estimate that ~ 10 — 15% of total events seen toward the M31 bulge area are 
caused by the Galactic halo MACHOs, assuming an all-MACHO halo. The first and 
second year of data observed by the MACHO group toward the Large Magellanic Cloud 
(LMC) gives tentative support to the hypothesis that MACHOs constitute a significant 
fraction of the Galactic halo (D. Bennett 1995, private communication). We test various 
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observables for their ability to distinguish the Galactic halo events from M31 events. 
We find that most high maximum-flux Galactic halo events can be isolated through a 
satellite-based measurement of the flux difference caused by the parallax effect. This 
observational strategy is discussed in § 4.1. Proper motion can be measured for ~ 5% of 
high signal-to-noise ratio (S/N) M31 events and can also be used to separate Galactic 
halo/M31 events. In addition, 1% of Galactic events can be isolated by determining the 
lower limit of the proper motion, fi\ ow . Hence, only a small fraction of events can be 
separated by using the proper motion and proper motion-limit measurements. 



2. The Fraction of Galactic Halo Events 



2.1. Optical Depth Estimate 



The mean optical depth to an M31 star lensed by the Galactic halo MACHOs is 

7"MW,halo = — o- / dD \PMW,halo(D \)D, (2.1.1) 

c Jo 

where c?M3i = 770 kpc is the distance to M31, and D = D Q \D\ S / D os , and D os , D \, and D\ s 
are the distances between the observer, source, and the lens. For the distribution of matter 
in the Galactic halo, we adopt a modified isothermal sphere with a core radius r c ^w; 

o^w Po( r o,MW + #o)/(<mw + r2 )i r < 200 k P c (o I 9) 

mw ' halo(rj " l0, r> 200 kpc, (2 - L2j 

where r is the distance measured from the Galactic center to a MACHO, the normalization 
is po = 7.9 x 10~ 3 M Q pc~ 3 , the core radius is r C) MW = 2 kpc (Bahcall, Schmidt, & Soneira 
1983), and Rq = 8 kpc is the solar galactocentric distance. The Galactic halo model is the 
same as that adoted by Griest (1991) for his computation of the optical depth toward the 
LMC except that the halo extends out to r = 200 kpc, and thus results in higher TMW,haio- 
The value of r C) MW is quite uncertain, but it hardly affects the optical depth determination 
because M31 is located ~ 119° away from the Galactic center. From equations (2.1.1) 
and (2.1.2), one finds the optical depth contribution by the Galactic halo MACHOs to be 

TMW,halo = 4.4 X 10~ 7 . 

For the M31 halo, we adopt a standard isothermal model with a core radius of the 

form 

PM31 ' hal ° " 4^ + r c V 31 ) (2 - 0) 
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because it is known that the rotation curve is flat at large radii up to r = 30 kpc from the 
M31 center with a rotation velocity of u c ,M3i = 240 km s _1 (Kent 1989a). The core radius 
is determined from the relation 

1)^ POO ( \_ '- \ 

M c ,M3i = M M3 i,disk + M M3 i,buige = J ^¥ 1 / drr 2 r — - 2 , (2.1.4) 

^ JO \r r +r cM31 j 

resulting in r Cj M3i = (2/-7r)(GM Ci M3i/^c M3i) = 6-5 kpc w hh the total mass of the 
luminous components taken to be M C) M3i = 13.2 x 10 10 M , where the total masses of 
the M31 bulge and disk are determined below. With this density model, the optical depth 
contributed by M31 halo lenses is computed similarly to equation (2.1.1), resulting in 
TM31,halo = 1-91 x 10~ 6 . 

On the other hand, the optical depth along a given line of sight caused by M31 
disk+bulge events is computed by 



AttG f d2 f D ° l \ f d ' 2 

TM31 = — 2~ / dD 0S pM3l{D s) / dD ol p M 3i {D Q \)D I dD os pM3l{D { 



(2.1.5) 



because theses populations work both as lenses and sources (self-lensing) . Here 
/>M3i = PM3i,disk + PM3i,buige, PM3i,disk and /9 M3 i,buige are the M31 disk and bulge densities, 
and d\ = dyizi — 30 kpc and d<i = dM3i + 30 kpc are the lower and upper boundaries of 
mass distribution along the line of sight. We model the M31 bulge as an oblate spheroid 
with its (unnormalized) matter density as a function of semimajor axis provided by 
Kent (1989b). The best-fitting value of axis ratio, c/a = 0.75, is found by fitting the 
computed surface brightness (including extinction) to a y-band image of M31 bulge 
(Han 1996). In the central region where the M31 bulge mass dominates, the rotation 
velocity is ~ 275 km s -1 . We thus normalize the total M31 bulge mass within r = 4 
kpc by M M 3i,bul g c = (275 km s _1 /«c,Mw) 2 ^MW,bulgc, where w c ,mw = 220 km s _1 and 
Mmw, bulge = 2 x 10 10 Mq are the adopted rotation speed of the Galaxy and Galactic bulge 
mass within r = 4 kpc (Zhao, Spergel, & Rich 1995). This normalization procedure results 

in -Mm 3 1, bulge 

= 4vr J °° drr 2 pM3i,buigc = 4.9 x 10 10 M , which is in a reasonable agreement 
with the determined value of 4.0 x 10 10 M by Kent (1989b). The M31 bulge is cut off at 
r ~ 8 kpc. The M31 disk is modeled by the double exponential disk, i.e., 

PM3i,d isk (^) = ^exp (-£) exp , (2.1.6) 

where So = 280 M & pc~ 2 is the normalization, and the radial and vertical scale 
heights are h z = 400 pc and = 6.4 kpc, respectively (Gould 1994c). This results in 
M M3 i,disk(# < 30 kpc) = 8.3 x 10 10 M . 

The resulting total optical depth distribution caused by M31 population is shown as 
a contour map in Figure 1 for various lines of sights. The projected position (x',y') is 
defined relative to the major and minor axes of the M31 bulge isophotes. We note that 
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the determined optical depth is subject to many uncertainties. The main uncertainty 
comes from the density distribution of lenses and sources. For example, Braun (1991) 
has argued that there is no need to invoke a massive dark halo component. Instead of a 
halo, his model has higher disk and bulge densities, c.f., Mm3i, bulge = 7.8 x 10 10 M and 
Mm3i disk = 12.2 x 10 10 Mq within 30 kpc from the center. Pixel lensing experiments 
may resolve the conflicts in the matter density distribution in M31 from the comparison 
between the theoretical predictions and future observational result. 



2.2. Event Rate Estimate 



Now we estimate the fraction of events caused by MACHOs in the Galactic halo when 
a pixel lensing search is carried out toward the M31 bulge from the ground. The event rate 
per angular area dfl is determined by 

^ = v^f^ E r o/W^Mi^); r = -u,T, (2.2.1) 

where 4>{Fq) is the luminosity function (LF) of the source stars (see § 3) and £ is the 
surface brightness in the M31 bulge area (see Han 1996 for detailed formalism and the 
surface brightness distribution). Here the maximum impact parameter, (3 m a.x(Fo), within 
which a lensing event with a source star of luminosity Fq can be detected for a certain 
required signal to noise ratio, (S/N) m i n , is obtained from the equation 

C[/WF ° )] - f^h^\ (S/N)l mV , (2.2.2) 



Pmax V 7rF 2 



a. 



where a is the mean photon detection rate averaged over the duty cycle of the telescope, 
*cxp = 1 day (see below), u) = (tc)" 1 is the inverse value of the median Einstein time scale, 
and Q ps f ~ tt^soc/ l n 4 is the angular area of the point spread function (PSF) of the source 
star with a seeing S ee- The values of C(P) sxg given in Figure 1 of Gould (1996). The 
signal suffers from additional noise due to the sky brightness and the time variation of the 
PSF by a factor 

2 2,2 

V = rj s ky + ri ps{ ; 

, ly = £±^, , ? , [ = !|^e(^)\ (2 ., 3) 

where the sky brightness is taken to be S s w ~ 19.5 mag arcsec -2 in / band, A6 see is the 
seeing difference between an image and the reference frame, and. 6?fi/ 

is the rate of photon 
detection at the fluctuation magitude of M31, i.e., fhi = 23.2 (Tonry 1991). 

The event rate, T, depends on the observational strategy. Here the event rate is 
determined under the following observational conditions. Event rates for various other 



-6- 



strategies are discussed in detail by Han (1996). We assume the experiment employs a 
2K x 2K CCD array with a pixel size of 0."25 on a 1-m telescope. The observation is 
carried out in / band, and the camera can detect 12 photons s _1 for an / = 20 star. The 
observation is carried out for an average 4 hrs per night during the M31 season, which is 
~ 1/3 of a year. The fractional seeing difference is kept to (A6 see /9 see ) = 5% and each 
image is obtained by subdividing the exposure time to prevent saturation of pixels. In the 
computation, the sub-integration time is ~ 20 min, and thus n^iv = 4 hr/20 min = 12. The 
images taken every day are combined together, and thus a = e x (4 hr/24 hr) photons hr _1 , 
where e is the photon detection rate of the camera. 

For the computation of event rate, it is required to compute the median time scale 
(t c ) = 1/tD [see equation (2.2.2)]. The value of (t e ) depends on the geometry, the mass 
function, and the transverse velocity of the lens system. The velocity distribution is 
modeled by a Gaussian with one-dismensional dispersion, a. The Galactic MACHO 
velocity distribution has o"MW,halo = 220/\/2 km s^ 1 . To account for the solar motion 
around the Galactic center, we increase the dispersion by 0Mw,haio = 250/\/2 km s _1 . 
By incresing the Galactic halo lens speed, one can simplfy the observer's motion to be 0. 
The M31 halo lens has 0"M3i,haio = ^c.msi/v 7 ^ ~ 170 km s _1 , while the disk and bulge are 
modeled by ctm31 = 156 km s^ 1 , based on the measured value by Lawrei (1983). With 
these adopted values of velocity dispersions the transverse velocity is determined by 



V = V] 



, Ad\ , ( As 
v s — + v, 



(2.2.4) 



where v s , vi, and v Q are the velocity of the source, lens, and observer, respectively. For 
Galactic halo events, one can approximate D \/D os ~ and D\ s /D os ~ 1. Similarily, 
D i/D os ~ 1 and D\ s /D os ~ for M31 halo and M31 disk+bulge self-lensing events. The 
values of a of the transverse speed distributions for events of individual populations are 
listed in Table 1. The resultant speed distributions of both M31 halo and disk+bulge 
events are similar, so we use a common distribution of <7M3i = 225 km s _1 . 

With the transverse speed distribution models, the time scale distribution of self 
lensing events (M31 disk+bulge events) for a constant mass, e.g., 1 M Q , is then obtained 
by 

/(f e ) = t dD osP (D os ) ( D ° l dDMDoi)D 1 ' 2 r dv5 (t c - y™^R \ vf{vh (2.2.5) 
Jdi Jd! JO \ CV J 

where S is the Dirac delta function. For the Galactic and M31 halo events, the time scale 
distribution can be approximated by 

/(O = /" maX dDMDo^ 2 [°° dv5 (t e - y™^R \ vf{v)i (2.2.6) 
id min Jo \ CV J 

where d m \ n = and d max = 200 kpc for Galactic halo events while <i m i n = dyizi — 200 kpc 
and d max = G?M3i for M31 halo events. The additional factors D 1 / 2 and v are included to 
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weight events by their cross section r c oc D 1 / 2 and transverse speed. Then the true time 
scale distribution is obtained by convolving f(t' e ) with the lens mass function /m(M); 

/(t e ) = | dMM l l 2 f M (M) J™ dt'j(t' c )5(t c - M 1 IH' C ). (2.2.7) 

Here once again the factor M 1//2 is included to weight by r c oc M 1 / 2 . For the mass of 
Galactic and M31 halo MACHOs, we adopt M = 0.1 M = constant, which is the upper 
mass limit of hydrogen burning. Note that t e oc M 1 / 2 for other masses. However, the 
majority of M31 disk+bulge events are expected to be caused by low mass stars, and 
thus we adopt a mass spectrum of local Galactic stars obtained from HST data (Gould, 
Bahcall, & Flynn 1996) for M31 disk+bulge lenses. The mass spectrum is shown in 
Figure 2 of Han & Gould (1996). To compare the time scale distributions for events from 
individual populations, we present f(t e ) for events at the position (x',y') = (1 kpc, 1 kpc) 
as representative distributions in Figure 2. Note that f(t e ) for M31 disk+bulge events 
depends on the field of observation; shorter t e toward the center because the average 
separation between lenses and sources ((As)) decreases. Events from the three populations 
have very similar distributions with (t e ) ~ 20 days and so the time scale cannot be used as 
a tool for the separation of events into different populations. 

Finally, the total event rate is then found by 

§-dn, (2.2.8) 

where the total angular area of the CCD is Occd = (500" x 500"), which is equivalent 
to ~ (1.87 x 1.87) kpc 2 . The event rate for each population as a function of (S/N) m - m is 
shown in Figure 3. The M31 disk+bulge and M31 halo events have a similar event rates, 
while Galactic halo events are ~ 3 — 4 times less frequent than events from either of the 
M31 populations, i.e., Galactic halo events constitute ~ 10 — 15% of the total. 

The fraction of Galactic halo events will increase by choosing a filed that is located 
away from the M31 bulge. As the field is further away from the M31 center, the optical 
depth for the Galactic halo events is virtually the same, but the source star number, iV*, 
and thus the total number of events, N tot , significantly decreases because 7V tot oc N*t. For 
example, the number of source stars in the field at (1 kpc, 1 kpc) is a factor ~ 45 larger 
than that of the field at (5 kpc, 5 kpc). Therefore, it would not be a good idea choosing a 
field too far away from the M31 bulge just to increase halo/M31 event ratio. 



3. Monte Carlo Simulation 



Now we return to the main concern of the article; what kind of information can one 
obtain from observables of individual events? To answer this question, we carry out a 
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Monte Carlo simulation of lensing events toward the M31 bulge. The following are the 
distributions from which lensing parameters are drawn; 



g M (M) = J M f(M')M'^ 2 dM', 
g v (v) = m(v')v'dv', 



9r c ( r e) 



J dD xp(D o[ )D 1 / 2 5 



1/2 



halo , (3 ^ 

1/21 



/ dD os p(D os ) J dDMD^D 1 /^ r c - (^^) ' j , M31, 

l > ^(/9) = [0,Anax], 

where [0, /3 max ] represents a uniform distribution in the range < (3 < (3 max . 

The source star LF is based primarily on the LF of Galactic bulge stars determined 
by J. Frogel (1995, private communication). However, this determination is incomplete 
at the faint end (Mj > 3.6). If it were adopted without correction, the fraction of 
luminous stars would be overestimated. The net effect would then be the overestimation 
of the event rate because events are more likely to be detected for luminous stars. We 
model the faint part of the LF by adopting the LF of stars in the solar neighborhood 
determined by Wielen, Jahreiss, & Kriiger (1983) for 3.6 < Mj < 8 and by Gould, 
Bahcall, & Flynn (1996) for 8 < Mj < 14. The corrected model LF yields a fluctuation 
magnitude of J2i <t>( F o,i) F o,i/ £i (f>( F o,i) F o,i = 23.5, which matches well with Tonry's (1991) 
determination of fhj = 24.3 before the extinction correction. The fraction of light from 
the stars in the corrected part of the LF is 14.7%, which is not negligible. The LFs before 
(thin line) and after the model correction (thick line) are shown in the first panel of Figure 
5. For each luminosity, the value of /3 max is obtained from equation (2.2.2). 

Detectable events are those that satisfy the condition (3 < /3 max . In addition, events 
must satisfy the following two conditions. The first one is that the effective duration of 
an event, t e g = (3t e , must be long enough for the reasonable construction of light curve, 
i.e., t c g > 0.5 day. Second, the angular size of the source star must be smaller than the 
angullar Einstein ring size, i.e., 9*/9 e < /3 max , to avoid the maximum flux suppression due 
to the finite size source effect (Nemiroff & Wickramashinghe 1994; Witt &; Mao 1994; 
Gould 1994a; Simmons, Newsam, & Willis 1995). Both Galactic and M31 halo events have 
very similar physical parameters. They have same masses (0.1 Mq) and similar transverse 
speeds (177 km s _1 versus 225 km s™ 1 ). In addition, due to the equivalence between D Q \ 
for Galactic lenses to D\ s of M31 halo events, vice versa, both population events have very 
similar time scales, as already shown in Figure 2, and Einstein ring sizes, r c oc y/MD/v. 



4. Separation of Halo/M31 Events 



Separating Galactic halo/M31 events is important for the better understanding of the 
global properties of the Galactic halo such as its flattening. Because Galactic halo events 
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account for only ~ 10 — 15% of the total, it is possible to use observations toward M31 
to constrain the Galactic halo only if the halo events can be identified. Sackett & Gould 
(1993) have proposed that one can constrain the shape of the Galactic halo by studying 
the optical depths toward LMC and Small Magellanic Cloud (SMC). To understand the 
role of M31 in determining the shape of the halo, we write the Galactic coordinates of the 
SMC, LMC, and M31 as (l,b) = (-57°, -44°), (-80°, -33°), and (121°, -21°). Thus, the 
lines of sight toward SMC, LMC, and M31 are progressively closer to the Galactic plane 
and further from the Galactic center. This implies that by adding the line of sight to M31, 
one increases the leverage on the shape of the Galactic halo relative to observing the LMC 
and SMC alone. 

Additionally, halo events are contaminants for the study of MACHOs in M31. If 
the Galactic halo is flattened (Dubinski & Calberg 1991; Katz 1991; Sackett & Sparke 
1990), the contamination of events toward M31 by the Galactic halo events would be more 
significant than our estimate. If one can separate a large fraction of M31 events from the 
Galactic halo events, a better study about the nature of M31 MACHOs will be possible. 
In this section, we investigate the methods and required observables that can be used to 
isolate the two populations. 

We begin with a brief discussion of the general observables in pixel lensing. Near the 
peak of a light curve, the flux is approximated by 

-Tmax = -T, feff = P*e, (4.1 J 



[l + ^Aeff) 2 ] 1 / 2 ' ^ (3 

since A — > 1/x and x = [(3 2 + {t/t c )\ l l 2 . Here -F max is the flux at the maximum 
amplification. The values of -F max and i e ff are directly measurable for individual moderately 
high amplification events. The distributions of t c g and F max are determined for events 
with (S/N) m i n = 20 and (S/N) m { n = 70 by the Monte Carlo simulation and are shown in 
Figure 4. Since we already know that Galactic and M31 halo events have similar physical 
paramaters, and so similar for F max and t e g, we show only distributons of Galactic halo 
and M31 disk+bulge events. We present their distributions at the representative positions 
of the M31 bulge area at (x', y') = (0.5 kpc, 0.5 kpc) {thick line) and (1.5 kpc, 1.5 kpc) {thin 
line). The distributions for the Galactic (and M31) halo and M31 disk+bulge self-lensing 
events are similar to each other, and thus neither quantity can be directly used for the 
separation of individual populations. Although t e can be measured for some of the high 
S/N (> 70) events, it is not in general a direct observable. However, one can infer t e from 
i e ff and F max provided that the color of the source star is known. The color of a source star 
can still be measured from the maximum flux since F max j/F max y = Fo,i/Fq,v and the 
event is achromatic. The approximate value of the luminosity, L, can then be determined 
using a color-mag relation. However, there can be different types of stars with the same 
color, e.g., lower main sequence stars and giants. Therefore, there will be a two-fold 
degeneracy in determining L for low S/N events. For high S/N events, t e (and so Fq) 
can be approximately determined from the light curve, thereby breaking the degeneracy 
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and allowing unambiguous determination of L from the color. Fortunately, the stellar 
population expected to be detected at low S/N are nearly all main sequence stars, so it 
should be possible to determine t e in nearly all cases. In Figure 5, we plot the LF of source 
stars (before amplification) for the expected events with different values of threshold S/N; 
(S/N) m i n = 10 and 70 in the second and third panels, respectively. At lower (S/N) m [ n , 
most detectable events are from main-sequence stars. As the (S/N) m i n increases, a 
significant fraction of events comes from upper main-sequence and clump giant stars. As 
shown in Figure 2, however, t e is not useful for separating halo/M31 events, either. Two 
other quantities can be measured under special circumstances: the displacement of the 
source star in the Einstein ring due to parallax effect and the proper motion of the lens, fx. 



4.1. Parallax Measurement 

One of the tools for the separation of both Galactic halo/M31 events is provided 
by the parallax effect. The principle is simple. By observing an event at two different 
positions, one from the ground and the other from a heliocentric satellite, the light 
curves display a displacement in the Einstein ring, Ax (Refsdal 1966; Gould 1994b). 
Then the flux measured from the ground will be different from that measured from the 
satellite for the Galactic halo events because the ground-based photometry is carried 
out when the event is near the maximum amplification, while the flux measured from 
the satellite, -F sa t, is not. However, the flux difference due to parallax effect will be 
extremely small for M31 self-lensing events due to their negligible amount of parallax 
[Ax ~ (D] B /D os )(AJJ /r e ) ~ 10~ 2 ]. Therefore, when the difference between the fluxes, 
AF, from ground and satellite-based photometry is detected, the event is almost certainly 
caused by a Galactic halo MACHO. It will be difficult to precisely measure Ax; the 
signal-to-noise ratio will be low and a typical event lasts only {t e s) ~ 1 day, which is too 
short for the construction of detailed light curve. However, measurement of AF will still 
be feasible, and this by itself allows the Galactic events to be recognized. 

Now, we answer the question: what fraction of Galactic halo events can be 
distinguished in this way? For this determination, we assume that the satellite photometry 
will be carried out with a camera on a 0.5 m telescope mounted on a heliocentric satellite. 
The CCD camera is exposed for At = 20 min and can detect 3 photons s _1 for an I = 20 
mag star. We assume that the seeing is slightly larger than the diffraction limit: 6 SCC = 0".5. 
The background surface brightness varies in the range £/ ~ 19 — 22 mag arcsec~ 2 in 
most of the M31 bulge area except for the part very close to the central region. The 
read-out noise is assumed to be negligible. Then the signal to noise ratio from the satellite 
observation is computed by 



(4.1.1) 
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where F sat is the measured flux from the satellite, the background flux is 
^back = F Ifi W-° ASi n psi , and F Ij0 = 3 x 10 8 photons s" 1 is the flux for 1 = 0. 

The flux difference due to parallax is computed as follows. The Earth-satellite 
separation vector projected onto the lens plane is Ar = R — (R • s)s, where R is the 
3-dimensional Earth-satellite vector and s is the unit vector toward M31 from the Sun. In 
scalar form, the separation is given by 

Ar = R(l -cos 2 /?c C cos 2 V) 1/2 , (4.1.2) 

where (3 ec = 33°.3 is the ecliptic latitude of M31 and -0 is the phase of the orbit (Gould 
1995). During an M31 season, ip varies in the range 60° < ip < 180°. The amount of 
displacement in Einstein ring due to the parallax is then calculated by 

Ar 

Sx = Ax cos <3? = — cos<3?; cos<I? = r-x, (4.1.3) 

r e 

where the projected Einstein radius is f c = r e (D os / D\ s ). For Galactic halo events, f e ~ r e 
since D os /D\ s — ► 1. The angle varies in the range [0, 2tt], i.e., random orientation, and 
the distribution of r c is obtained from equation (3.1). Then the flux difference is computed 
by 

AF = F max - F sat = [A(/3) - A((3 + 5x)} F . (4.1.4) 
In Table 2, we present the numbers, -/V se p, and fractions, A se p / /7V~]y['\y 1 haio> 

of Galactic 

halo events that can be separated from M31 events by measuring AF for various values 
of F max . Both values are determined under the criteria of (S/N)$ x > 3 and 5 (i.e., 3 
and 5 a levels). The numbers of events due to the Galactic halo, A^MW.haioj and due to 
all populations (M31 halo, disk+bulge, and Galactic halo), iVtot, are those that can be 
detected with S/N > 20 from ground observations. For event with -F max < 23 mag, nearly 
all Galactic halo events can be separated, while the fraction decreases significantly for low 
F max events because the stars are too faint to be seen. (To avoid confusion, we note that 
these faint events can still be detected from the ground because the image combination 
process yields a higher S/N.) Therefore, the best target for AF measurement would be 
high 

Fma,x events. Since the photometry can be done with a relatively short exposure time, 
At ~ 20 min, it can be carried out as a part of a larger project whose principal mission is 
parallax measurements of events toward the Galactic bulge and LMC. 



4.2. Proper Motion 

Proper motion of MACHOs can also be used for the separation of Galactic halo/M31 
events. For a classical lensing event in which a MACHO passes over or very close to 
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the face of a star, the part of the star closer to the lens is amplified more than other 
parts because the source star is not a perfect point source: differential amplification. By 
measuring the amount of deviation of the light curve from that of point source, one can 
measure the ratio 

= |, (4.2.1) 

where 6>* = R*/D os is the angular size of the source. Since 9* can be determined from 
the measured color (and Stefan's law), one can uniquely determine the 9 e and the proper 
motion of the MACHO, fi = 9 e /t e (Gould 1994a). For pixel lensing events in which the 
lens transits the source, it is also possible to measure the proper motion. However, the 
process is somewhat different from that in classical lensing. The source cannot be resolved, 
so it is impossible to determine 9* from Stefan's law. Instead, one would estimate R* from 
the observed color and the color-mag relation, and then compute 9*. Generally, it is also 
not possible to measure t e . However, from the deviation of light curve from its point source 
form, one can determine t*, which is the time it takes for the lens to cross the stellar 
radius. The proper motion is then given by ^ = 9*/t*. Once /i is measured, one can easily 
distinguish halo/M31 events because the values of \i are quite different for the Galactic 
halo and M31 events. However, the fraction of events for which proper motions can be 
measured is small. For M31 disk+bulge self-lensing events with S/N > 50, the fraction 
is ~ 5% and similar to M31 halo events. Since the fraction scales inversely with 9 e , the 
Galactic fraction is smaller by a factor ~ 10~ 2 , and is therefore completely negligible. 

In the more usual case when the lens does not transit the source star [j3 > (9*/9 c )], 
one can still set the lower limit of the proper motion of individual events by 

Q 

A* > Wow = t^- (4.2.2) 

teft 

Then a Galactic halo event can be isolated whenever ni ow is too high to be consistent with 
M31 self-lensing. We determine the fraction of Galactic halo events that can be isolated 
in this way by setting the upper limit of the M31 lens population at 1500 km s -1 <f^ 31 . 
Unfortunately, only ~ 1% of Galactic halo events can be isolated using proper motion 
lower limits. 



This work was supported by a grant AST 94-20746 from the NSF. 



-13- 



REFERENCES 

Bahcall, J. A., Schmidt, M., k Soneira, R. M. 1983, ApJ, 265, 730 

Baillon, P., Bouquet, A., Giraud-Heraud, Y., k Kaplan, J. 1993, A&A, 277, 1 

Braun, R. 1991, ApJ, 372, 54 

Ciardullo, R., Tamblyn, P., k Phillip, A. C. 1990, PASP, 102, 1113 

Crotts, Ap. P. S. 1992, ApJ, 399, L43 

Dubinski, J., k Calberg, R. G. 1991, ApJ, 342, 212 

Gould, A. 1994a, ApJ, 421, L71 

Gould, A. 1994b, ApJ, 421, L75 

Gould, A. 1994c, ApJ, 435, 573 

Gould, A. 1995, ApJ, 441, L21 

Gould, A. 1996, ApJ, 470, 000 

Gould, A., Bahcall, J. A., k Flynn, C, 1996, ApJ, 465, 000 

Griest, K. 1991, ApJ, 366, 412 

Han, C, k Gould, A. 1996, ApJ, 467, 000 

Han, C. 1996, ApJ, submitted 

Katz, N. 1991, ApJ, 368, 325 

Kent, S. M. 1989a , PASP, 101, 489 

Kent, S. M. 1989b, AJ, 97, 1614 

Lawrie, D. G. 1983, ApJ, 273, 562 

Nemiroff, R. J., k Wickramasinghe, W. A. D. T. 1994, ApJ, 424, L21 

Refsdal, S. 1966, MNRAS, 132, 101 

Sackett, P. D., k Sparke, L. S. 1991, ApJ, 371, 443 

Sackett, P. D., k Gould, A. 1993, ApJ, 419, 648 

Simmons, J. F. L., Newsam, A. M., k Willis, J. P. 1995, MNRAS, 276, 182 
Tomaney, A. B., & Crotts, A. P. S. 1996, ApJ, submitted 
Tonry, J. L. 1991, ApJ, 373, LI 

Wielen, R., Jahreiss, H., k Kriiger, R. 1983, IAU coll. 76: Nearby Stars and the Stellar 
Luminosity Function, eds. A. G. D. Philip k A. R. Upgren, 163 

Witt, H. J., k Mao, S. 1994, ApJ, 430, 505 

Zhao, H., Spergel, D. N., k Rich, R. M. 1995, ApJ, 440, L13 



This preprint was prepared with the AAS IAT^X macros v3.0. 



-14- 



FIGURE CAPTION 

Figure 1: The optical depth distribution toward the M31 bulge. The contours are 
drawn in units of 10~ 6 . In the computation, the uniform contribution by the M31 halo, 
TM3i,haio = 1-91 x 10~ 6 , is included. The total optical depth both by M31 and Galactic 
halo lenses is obtained by adding Thaio to the values marked on the contour, where 
Thalo = 4.4 x 1CT 7 . With the adopted distance to M31 of cZm31 = 770 kpc, 1 arcsec 
corresponds to ~ 3.73 pc. 

Figure 2: The Einstein time scale distributions, f(t c ), for the Galactic halo, M31 
self-lensing, and M31 halo events at the position (x',y') = (1 kpc, 1 kpc). The time scale 
distributions for individual populations are very similar one another. 

Figure 3: The event rates, T, for individual population events as a function of the 
threshold signal-to-noise ratio, (S/N) m i n . The observational conditions are described in § 
2.2. 

Figure 4: The effective time scale and maximum flux distributions for (S/N) m [ n = 20 and 
70 for the Galactic halo (dotted lines) and M31 disk+bulge (solid lines) events. Because 
the distributions are location dependent, those at (x',y r ) = (0.5 kpc, 0.5 kpc) (thick lines) 
and (1.5 kpc, 1.5 kpc) (thin lines) are presented as representative distributions. 

Figure 5: The LFs of detectable events with different values of (S/N) m \ n . The original 
LF, modeled by combining the Galactic bulge LF (J. Frogel, private communication), thin 
line, and the LFs of Gould et al. (1996) and Wielen et al. (1983), thick line, is shown in 
the first panel. The (unamplified) LF of detectable events with (S/N) m i n = 10, and 70 are 
shown in the second and third panels. 
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TABLE 1 
TRANSVERSE VELOCITY 



event 


dispersion 


(- 2 ) 


population 


(km s _1 


) 2 


Galactic halo events 


177 2 + 2 = 


= 177 2 


M31 halo events 


170 2 + 156 2 = 


= 230 2 


M31 (disk+bulge) events 


156 2 + 156 6 = 


= 221 2 



NOTE. — The transverse speed distribution for each population event. Each speed has a Gaussian 
distriburion with a dispersion listed. 
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TABLE 2 

GALACTIC HALO EVENTS SEPARABLE BY PARALLAX EFFECT 



events 


N tot 


-^MW.halo 




(S/N) Sx > 3 




{S/N) Sx > 5 


with F max 






iVsep 


-^Vsep/-^MW,halo 


iVsep 


^Vscp/^MW.halo 


< 21 mag 


2.62 


0.13 


0.13 


100% 


0.13 


100% 


< 22 mag 


7.42 


0.37 


0.36 


98% 


0.36 


97% 


< 23 mag 


18.95 


0.96 


0.92 


96% 


0.84 


87% 


< 24 mag 


41.59 


2.38 


1.36 


57% 


0.86 


36% 


< 25 mag 


65.31 


5.19 


1.40 


27% 


0.94 


18% 


total 


95.00 


10.0 


1.70 


17% 


1.20 


12% 



NOTE. — The numbers, N scp , and fractions, N sep /N M w,halo, of Galactic halo events that can be 
separated from M31 events by measuring AF for various values of F miH . Both values are determined under 
the criteria of (S/N)s x > 3 and 5 (i.e., 3 and 5 a levels). The numbers of events due to the Galactic halo, 
A^MW.haio, and due to all populations (M31 halo, disk+bulge, and Galactic halo), JVtot, are those that can 
be detected with S/N > 20 from ground observations. For event with F max < 23 mag, nearly all Galactic 
halo events can be separated, while the fraction decreases significantly for low F max events. 
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